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We present a theory to study gas molecules adsorption on armchair graphene nanoribbons (AG- 
NRs) by applying the results of ab initio calculations to the single-band tight-binding approximation. 
In addition, the effect of edge states on the electronic properties of AGNR is included in the cal- 
culations. Under the assumption that the gas molecules adsorb on the ribbon sites with uniform 
probability distribution, the applicability of the method is examined for finite concentrations of ad- 
sorption of several simple gas molecules (CO, NO, CO2, NH3) on 10-AGNR. We show that the states 
contributed by the adsorbed CO and NO molecules are quite localized near the center of original 
band gap and suggest that the charge transport in such systems cannot be enhanced considerably, 
while CO2 and NH3 molecules adsorption acts as acceptor and donor, respectively. The results of 
this theory at low gas concentration are in good agreement with those obtained by density-functional 
theory calculations. 



I. INTRODUCTION 

In the past two decades carbon nanostructures, such as 
carbon nanotubes (CNTs) and graphene, have attracted 
much attention due to their novel fundamental properties 
and possible applications in future carbon-based nano- 
electronics [l[ . An important feature of these materials is 
their sensing property and on this basis, gas sensors 
pH sensors [j], and biosensors [H] have been proposed. 
Among various types of sensors, most of the experimental 
ill and theoretical [l4-22] studies have been focused 
on the gas sensors characteristics of the carbon nanos- 
tructures due to their promising applications in various 
fields such as electronics, agriculture and medicine 

It has been demonstrated that the CNTs can be used 
as chemical sensors for detecting very small concentra- 
tion of NO2, NH3 and other gases with high sensitivity 
at room temperature The interaction of graphene 

with chemical environments and the possibility of detect- 
ing individual molecules have also been reported fTol - fl~2l j . 
The results indicate that the sensing mechanism is based 
on changes in the charge carrier concentration induced 
by gas molecules adsorbed on the graphene surface and 
acting as donors or acceptors. For instance, Leenaerts et 
al [2l| . based on the first-principle calculations, showed 
that H 2 and NO2 behave as acceptors while NH 3 , CO 
and NO act as donors. They found that, molecular dop- 
ing, i.e. charge transfer between the adsorbates and the 
graphene surface, is almost independent of the adsorption 
site, which can be attributed to the translational symme- 
try of the system. In fact, due to the flat structure and 
thus larger accessible surface area, graphene as well as 
graphene nanoribbons (GNRs) may act better than any 
other carbon-based materials as gas sensors. 

The special electronic behavior of GNRs corresponds 
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to the two typical topological shapes of the carbon atoms 
on their edges, namely armchair and zigzag. These rib- 
bons can be either semiconducting with a size dependent 
gap or metallic 23] . The first-principles calculations and 
the results of simple tight-binding approximation show 
that all AGNRs with edge deformation are semiconduct- 
ing with a finite band gap 0, HH . Due to the semicon- 
ducting feature of the AGNRs, their sensing properties 
were investigated recently, based on the tight-binding ap- 
proximation [2^, [53] and the ab initio calculations [28j . 
By changing the values of hopping integrals and atomic 
on-site energies for simulating the gas adsorption on the 
edges of an AGNR and using the coherent potential ap- 
proximation for studying the effect of finite concentration 
of gas molecules, the sensing properties of the AGNR 
were discussed in Ref. Based on density functional 

theory (DFT) calculations, Huang and co-workers (28j 
studied the adsorption of CO, NO, N0 2 , 2 , N 2 , C0 2 , 
and NH3 on the edges of an AGNR in the presence of 
dangling bond defects. The gas adsorption only around 
the defect sites, centered at the edges, was considered. 
They found that the adsorption of C0 2 and 2 molecules 
changes the AGNR to p-type semiconductors, while NH3 
adsorption changes the system to an n-type semiconduc- 
tor. 

We should note that although the ab initio techniques 
based on DFT obtain a good degree of accuracy to de- 
scribe the electronic properties of pristine materials, it 
is somewhat restricted when we consider disordered sys- 
tems such as gas sensors. The reason is that, in a real dis- 
order system such as gas adsorption on GNR, one should 
consider all of the configurations and obtain a physical 
quantity by averaging over all possible adsorption config- 
urations, while in the DFT calculations only one config- 
uration is usually considered. 

In order to fully understand the effects of adsorbed 
molecules on the electronic properties of AGNRs, we 
should consider all possible configurations of gas adsorp- 
tion, i.e. the variation of the position of the gas molecule 
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FIG. 1: (Color online) (a) Schematic view of the pristine 
10-AGNR with modification in bond lengths at the edges 
(blue lines). The arrow shows the periodic direction and the 
red (black) dotted rectangle indicates an unit cell (a super- 
cell). (b)-(d) The structure around adsorbed molecules after 
the adsorption of diatomic, triatomic and quadratomic gas 
molecules, respectively. Only a single atom (red atom) of 
each molecule is connected to the ribbon. 



across the ribbon width. Therefore, using the values of 
bond lengths between the AGNR and the molecules, and 
also the values of inter-atomic distances obtained from 
the DFT calculations [28[ and using a scaling rule of the 
tight-binding hopping integrals, we study the effects of 
CO, NO, CO2 and NH 3 adsorption on the average local 
density of states (LDOS) at different gas concentrations 
within the single-band tight-binding approximation. The 
paper is organized as follows. In Sec. II, we present the 
model and formalism to obtain the Green's function and 
hence the average LDOS of the system. The total Green's 
function is obtained by means of the iterative procedure 
[29} , which provides a quick way of evaluating the surface 
and bulk Green's functions. Numerical results and dis- 
cussions for the electronic properties of AGNR with gas 
adsorption are presented in Sec. III. We conclude our 
findings in Sec. IV. 



ribbon consists of N a = 10 dimer lines across the ribbon 
width, as shown in Fig. 1(a), and following conventional 
notation [3(J, the ribbon is referred as 10-AGNR. We as- 
sume that the dangling er bonds at the edges have been 
passivated by hydrogen atoms (not shown), hence we ex- 
pect the bond lengths between carbon atoms at the edges 
differ from those in the interior of the ribbon (see the blue 
bonds in Fig. 1(a)). In addition, the gas molecules can 
be diatomic, triatomic or quadratomic and we assume 
that each gas molecule adsorbs to only one carbon atom 
(i.e., single contact) as shown in Fig. l(b)-l(d). In such 
a device, the electronic structure of entire system can 
be described by a single-electron Hamiltonian in a basis 
of localized atomic orbitals for the atoms in the molecule 
and 7r-orbitals in the 10-AGNR. We use the tight-binding 
approximation for both the ribbon and the gas molecules. 
The Hamiltonian of the system can be expressed as 
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where 'Kr corresponds to the pristine 10-AGNR (before 
gas adsorption) and can be written as 
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where the operator c\ creates (annihilates) an elec- 
tron at site i of the ribbon, ti is the on-site energy and 
tij is the hopping integral. Since there is no spin effect, 
we have suppressed the spin indices here. 

Because of the hydrogen passivation of the edge carbon 
atoms, both ii and ty of the edges can change. Here, 
we only consider the variation in the hopping integrals 
and we set = as an origin of the energy. It is re- 
ported that the hopping integral between 7r-orbitals at 
the edges of AGNRs increases about Siij = 12%tn [24j . 
Therefore, we set £y = 1.12 for the edges and Uj = ta 
with tR = 2.66 eV for the interior carbon atoms [3lj ]. 
Also, the Hamiltonian of the adsorbed molecules, "Km, 
is obtained by summation over the Hamiltonian of the 
isolated molecules (h# for molecule £) as: 



"Km = ^2 hi 



(3) 



with 



II. MODEL AND FORMALISM 

We consider an AGNR with gas molecules randomly 
adsorbed on its surface. The system is described as an 
infinite stack of supercells with nearest-neighbor interac- 
tion. This means that we transform the original system 
into a linear chain of supercells. Each supercell may con- 
tain several unit cells and the gas molecules adsorb to 
the AGNR sites with uniform probability. Under the as- 
sumption that only a single molecule adsorbs to each su- 
percell, the number of supercell carbon atoms, N s , deter- 
mines the gas molecule concentration as x = l/N s . The 



he = 



E 

(ni,mi) 
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where d' n {dm) creates (annihilates) an electron on 
atomic orbital n of the molecule £ and for simplicity, the 
atomic basis orbitals are taken to be orthonormal. Ex- 
tension to the non-orthogonal bases is straightforward in 
principle. The quantities e ne and t neme are the on-site en- 
ergy and hopping integral parameters for the atoms in the 
gas molecule. These parameters depend on the type and 
the geometrical arrangements of atoms in the molecule. 
For instance, CO2 molecule consists of two type atoms, 
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FIG. 2: (Color online) Electronic band structure of 10- 
AGNRs (a) before and (b)-(e) after gas molecule adsorption: 
(b) CO, (c) NO, (d) C0 2 , (e) NH 3 . In the case of adsorption, 
each unit cell contains a single molecule connected to one of 
the ribbon edges. F and X points correspond to the wave 
numbers k = and k = ±7r /\fia in the first Brillouin zone of 
the ribbon, respectively. 



carbon and oxygen, and there is no bonding between two 
oxygen atoms in the molecule, i.e. io-o = 0. Hence, we 
may have different values for the t ni and t neme . Finally, 
the interaction between the pristine 10-AGNR and the 
gas molecules is written as a sum over each molecule: 



with 
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where the indice ni labels the atomic site n of the 
molecule £ that is adsorbed on site i = it of the ribbon, 
and Un = t is the hopping integral between the ribbon 
and the molecule. N m is the number of atoms at each 
molecule. 

Now we can calculate the Green's function of the whole 
system to obtain the average LDOS. For this purpose, 
as mentioned above, we divide the system into an infi- 
nite series of supercells. The supercells are labeled by 
I = 0, ±1,±2, •■• and we focus on one of them, say 
£=0, and consider the supercells £ = — 1,— 2, ■•• and 
I = +1, +2, • • • as neighboring supercells at the left and 
right, respectively. Based on the iterative procedure 
introduced by Lopez Sancho et al. (29|, the Green's 
function of each supercell (in the presence of the gas 
molecules), due to the periodicity of the system, can be 
expressed as 

6(E) = (e - H 00 - H 01 f - H^T)- 1 , (7) 

where e = (E + iff) I with rj being a positive infinitesimal 
number and / is the identity matrix. Hqq represents the 



coupling between atoms within each supercell and Hqi 
stands for the interaction between two nearest-neighbor 
supercells. We note that, both are matrices of dimension 
(N m + N s ) x (N m + N s ). If we designate the Hamiltonian 
of a single supercell (say, £ — 0) of the pristine 10-AGNR 
by , then Hqq can be written in matrix form as 



and 
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where 6 is a (N m + N s - 10) x (JV m + N s - 10) zero ma- 
trix and O'Cft 1 ^ corresponds to the interaction between 
nearest-neighbor zigzag lines (perpendicular to the peri- 
odic direction) in the £ = and £ = 1 supercells. Also, 

T and T are the transfer matrices which can be easily 
computed using the matrices Hqq and Hqi through the 
iterative procedure [29l| . 
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where Pk and can be determined from the following 
recursion relations: 
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The process can be performed until the desired precision 
is achieved, that is, until P qi P q < 8 with 8 arbitrarily 
small. 

To study the electronic properties of the system, we ob- 
tain the diagonal matrix elements of the supercell Green's 
function at the ribbon sites and average over all config- 
urations. Therefore, the the average LDOS per carbon 
atom can be expressed as 



(D{E)) 



1 



(12) 



where Gi (i,i;E) corresponds to the diagonal matrix ele- 
ments of G(E) at the ribbon sites, when the gas molecule 
is contacted to the ribbon at site i — iy. 

III. RESULTS AND DISCUSSION 

We now use the model described above to study gas 
sensing properties of the system. The results of numer- 
ical computations will be presented for CO, NO, CO2, 
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FIG. 3: (Color online) Average LDOS of 10-AGNR with CO 
adsorption for two values of the gas concentration, in compar- 
ison with the LDOS of pure system. The dashed line shows 
the original ribbon Fermi energy. 



and NH 3 gas molecules adsorbed on the 10-AGNR. These 
simple molecules have been chosen due to their different 
electronic properties in the adsorption process on the rib- 
bon. In order to obtain the hopping parameters related 
to the adsorbates, we use of the scaling form of the tight- 
binding hopping integrals as [32j 



tr 
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where dn = 1.42 A is the bond length between carbon 
atoms in graphene , while d a p can be the bond length 
between the ribbon and the molecule, i.e. a £ ribbon and 
j3 G molecule, or between the atoms in the molecule, i.e. 
a and /3 £ molecule. We should note that, this formula 
is a valid approximation for the variation of the hopping 
parameter for a bond between two specific atoms when 
the distance between the atoms is slightly varied. In 
spite of this, we believe that, by scaling all the hopping 
parameters with this formula is able to qualitatively 
reproduce the electronic states of adsorbates similar to 
the ab initio results. To calculate the exact values for 
t a p, we need d a p parameters. Accordingly, we use the 
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FIG. 4: (Color online) The same as Fig. 3 but for NO ad- 
sorption. 



DFT results for these bond lengths to investigate the 
adsorption effects of the above-mentioned molecules [28[ ■ 
Moreover, the exact values for e n , are found with varying 
the on-site energies to obtain the desired LDOS which is 
consistent with the reported DFT results [2^|. Below 
we present the values of these parameters and discuss 
the adsorption effects of each adsorbate on the ribbon 
separately. Note that, we can also use the results of Ref. 

21 for comparison, in addition to Ref. (28{, because in 
this study the gas molecules can be adsorbed by each 
carbon atom across the ribbon width, contrary to Ref. 

281. 



A. CO on 10-AGNR 

We first emphasize that, according to the DFT cal- 
culations [ll|, the size of charge transfer between the 
molecules and the ribbon depends on the orientation 
of the molecules with respect to the ribbon surface. 
Moreover, the highest occupied molecular orbital in CO 
molecule is located on the C atom. Therefore, CO 
molecule is attached to a single atom of the ribbon via 
its carbon atom with the bond distance dc-c = 1-35 A, 
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while the bond length of the adsorbed CO is dc-o = 1-18 
A 28]. Inserting these values into Eq. [T3J we ob- 
tain r = 1.11 tn, and ic-o = 1.45 i^. Also, the ap- 
propriate on-site energies are given by ec = — 1.40 £r 
and eo = — 1.40 i^. Hence, the effect of this adsorbed 
molecule on the electronic structure of the system can be 
calculated numerically. 

The band structure and the average LDOS of 10- 
AGNR with CO adsorption are shown in Figs. 2(b) and 
3, respectively. The band structure, which corresponds 
to gas concentration x = 0.05, has been calculated for 
a configuration in which the molecule is connected to 
one of the ribbon edges. In comparison with Fig. 2(a) 
for the pristine 10-AGNR, it induces an impurity sate 
in the band gap, slightly above the Fermi energy of the 
clean system. To see the effects of all configurations in 
the presence of finite concentrations of gas molecule, we 
have plotted the average LDOS for two concentrations; 
x = 0.05 which indicates that each supercell with one 
unit cell contains a single CO molecule, as shown in Fig. 
3(a), and x = 0.0083 which indicates that each super- 
cell with six unit cells contains a single CO molecule as 
shown in Fig. 3(b). 



B. NO on 10-AGNR 

In adsorption of NO molecules on the 10-AGNR, only 
N atoms are attached to the ribbon [28|, with the bond 
distance <fc-N = 1-43 A, while the bond length of the 
adsorbed NO is g?n-o = 1-24 A. Accordingly, using Eq. 
[13] the hopping energies are given as r = 0.98 tR and 
iN-O = 1-31 tR. In addition, the appropriate values for 
the on-site energies are eN = —1.42 tn and eo = —1.35 tR. 
In this case, a single NO molecule on one of the ribbon 
edges induces an impurity sate in the band gap, slightly 
below the Fermi energy of the pure system [see Fig. 2(c)]. 
The results of LDOS clearly indicate that the values of 
charge density of the pristine 10-AGNR are modulated by 
averaging over all adsorption configurations. This type of 
adsorbate, similar to CO molecule, does not change the 
position of the bottom (top) of the conduction (valence) 
band. 

From the average LDOS in Figs. 3 and 4, we notice 
that the adsorption of both CO and NO molecules by 
the AGNR decreases the original band gap, especially 
at high x, and may affect the optical properties of the 
system. At low gas concentration, such as x — 0.0083, 
the density of states of gas molecules consists of a 6- 
function peak near the center of band gap, while at high 
x, such as x = 0.05, a narrow band appears instead of a 
level. However, this narrow band is far from the bottom 
(top) of the conduction (valence) band and one expects 
that such a gas adsorption does not enhance the charge 
transport of the system. 

In the case of CO and NO adsorption, we expect that 
the Fermi level does not change considerably due to the 
position of molecular level inside the band gap with re- 
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FIG. 5: (Color online) The same as Fig. 3 but for CO2 ad- 
sorption. 



spect to the edges of valence and conduction bands. As 
we know, in a non-degenerate n-type (p-type) semicon- 
ductor, the Fermi energy is rather far from the edge of 
the conduction (valence) band. This feature indicates 
that the 10-AGNR with adsorbed molecules like CO and 
NO behaves as a non-degenerate semiconductor. 



C. C0 2 on 10-AGNR 

The DFT results of CO2 adsorbed on the ribbon show 
that the molecule sits 1.51 A away from the ed ge c arbon 
atoms and the length of C-0 bonds is 1.26 A [28|]. Ac- 
cordingly, the calculated hopping energies are given as 
t = 0.89 £r and tc-o = 1-27 tR. Moreover, we choose 
ec = 1-75 tR and eo = 1-40 tR as the suitable values for 
the on-site energies. Using these parameters, the band 
structure for the configuration where the molecule sits 
on one of the edges is given by Fig. 2(d). An interesting 
feature of this figure is that the molecular states form an 
impurity band in the original band gap with width 0.3 
eV. We see that, due to the type of charge exchange be- 
tween the gas molecule and the ribbon, the energy levels 
of the molecule are close to the edge of the valence band. 



6 



0.30 




-10-AGNR with NH 3 adsorption 
- Pristine 10-AGNR 

x=0.05 



(a) 



0.00 





-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

Energy (eV) 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

Energy (eV) 

FIG. 6: (Color online) The same as Fig. 3 but for NH 3 ad- 
sorption. 



Also, the average LDOS analysis clearly show that, at 
x = 0.05 (Fig. 5(a)) the impurity band overlaps with 
the valence band and consequently the band gap reduces 
significantly. On the other hand, for low enough doping, 
such as x — 0.0083, the impurity states become more 
localized on the gas molecules, the adsorbate bandwidth 
decreases and a gap opens between impurity and the top 
of the valence bands [see the arrow in Fig. 5(b)]. This 
type of molecular doping can enhance the charge trans- 
port in this system and hence the CO2 molecule will act 
as an acceptor at finite x and temperature. 



D. NH 3 on 10-AGNR 

The ab initio studies of gas adsorption on graphene 
[2l| and 10-AGNR [28[ indicate that in the ammonia 
molecule, the N atom is closer to the surface than the 
H atoms. The bond lengths of the adsorbed NH3 are 
given as dc-N = 1.49 A and c?n-h = 1.03 A (Hi. 
Therefore, the hopping parameters are r = 0.9 tn and 
£n_h = 1.9 tfj, and the on-site energies are chosen to 
be cn = — 3A5tn and en = — 2.65 £r. Fig. 2(e) shows 
the electronic band structure of the system where the 



molecule sits on one of the ribbon edges. In this case the 
impurity band is close to the bottom of the conduction 
band with width 0.38 eV. 

We see in Fig. 6(a) that in the case of £=0.05, by aver- 
aging over all gas adsorption configurations, the impurity 
band overlaps with the conduction band and the band 
gap decreases considerably. In addition, at £=0.0083 as 
shown in Fig. 6(b), the impurity bandwidth decreases 
and at low enough x the impurity band is separated from 
the conduction band and a gap gradually opens [see the 
arrow in Fig. 6(b)]. This indicates that NH3 molecule 
acts as a donor, which is in good agreement with experi- 
mental [l(| and theoretical studies of NH3 adsorption on 
graphene surface [H| and 10-AGNR H|. 

From the above discussions we conclude that, C02 or 
NH3 molecules adsorption shifts the Fermi energy to- 
wards and even into the original valence or conduction 
band, and the system exhibits p-type or n-type semicon- 
ducting behavior. On the other hand, when the Fermi 
energy is close to or even higher (lower) than the bot- 
tom (top) of the conduction (valence) band, we have a so 
called degenerate semiconductor. Therefore, by looking 
at the Figs. 5 and 6, one can conclude that the adsorption 
of CO2 or NH3 molecules converts the electronic struc- 
ture of 10-AGNR as a degenerate semiconductor. We 
should note that, this feature which strongly depends on 
the gas concentration has not been investigated in the 
reported DFT calculations. 

By comparing Figs. 2(b)-2(e) which correspond to the 
gas adsorption configuration on one of the ribbon edges 
with Figs. 3-6 which correspond to the average of LDOS 
over all configurations, we conclude that the interaction 
between gas molecules and the carbon atoms at the rib- 
bon edges is much stronger than that of the interior car- 
bon atoms, which is consistent with the DFT results [28l . 



IV. CONCLUSION 

In this research, using the parameters obtained from 
DFT results, we have modeled the adsorption of gas 
molecules CO, NO, C0 2 , and NH 3 on the 10-AGNR. 
Under the assumption that the variation of bond lengths 
between carbon atoms is only localized on the edges of 
ribbon, and including this effect in our calculations, we 
studied the influence of finite concentrations of the ad- 
sorbates on the average LDOS. In addition, we assumed 
that only one molecule can occupy the ribbon sites and 
there is no interaction among the adsorbed molecules. 

The results show that, although all types of gas 
molecules can influence the electronic properties of the 
system, the average LDOS varies from one type to other. 
In the cases of CO and NO molecules adsorption, the im- 
purity states are very localized at the center of original 
band gap. Therefore, adsorption of these gas molecules 
on 10-AGNR cannot notably affect the charge transport 
through such systems. In the cases of CO2 and NH3 
adsorption, however, the system exhibits p-type and n- 
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type semiconducting behavior, respectively. Despite the supported by Payame Noor University grant, 
simplicity of the model, reasonable agreement with DFT 
results is demonstrated. The results can be useful in 
designing future gas sensors and p-type or n-type semi- 
conductors. 
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